Time-dependent Frohlich transformation approach for 
two-atom entanglement generated by successive passage through a cavity 
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Time-dependent Frohlich transformations can be used to derive an effective Hamiltonian for a 
class of quantum systems with time-dependent perturbations. We use such a transformation for a 
system with time-dependent atom-photon coupling induced by the classical motion of two atoms 
in an inhomogeneous electromagnetic field. We calculate the entanglement between the two atoms 
resulting from their motion through a cavity as a function of their initial position difference and 
velocity. 
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I. INTRODUCTION 

Canonical transformations have been widely used in 
condensed matter physics [l|, Q to derive effective Hamil- 
tonians by eliminating degrees of freedom with low- 
energy excitations [3]. One of the most well-known ap- 
plications was Frohlich's derivation Q of an attractive 
electron-electron interaction from the original electron- 
phonon interaction. 

Frohlich's approach (which was also studied by Naka- 
jima was to apply a unitary transformation H' = 
eyip(—gS)Hexp(gS), defined by an anti-Hermitian oper- 
ator S, to the Hamiltonian H = H Q + gH , where H° 
describes noninteracting electrons and phonons. His goal 
was to treat the electron-phonon interaction gH . The 
small parameter g is introduced to stress that the inter- 
action Hamiltonian is perturbative compared with the 
free Hamiltonian and can be set to 1 during/after the 
calculation. The transformation can be evaluated order 
by order, 



H' = H° + gH 1 + g[H°, S] + g 2 [H\S] 
+ 3l[[H ,S},S}+O(g 3 ). 



(1) 



Frohlich eliminated the term linear in g by requiring 

H 1 + [H°,S] = (2) 

to get the generator S. In the following, we will use the 
term Frohlich transformation for a canonical transfor- 
mation that fulfills Eq. @. The transformation leads to 
a phonon-induccd interaction among electrons with the 
potential V = g 2 ([H 1 , S}/2) by averaging over the low- 
energy phonon states. This effective Hamiltonian can be 
attractive or repulsive with a singularity at the energy 
shell. 

Since the interaction between atoms and electromag- 
netic field is similar to the electron-phonon coupling, it 
is natural to try to use Frohlich transformations in quan- 
tum optics [6]. In the large-detuning limit, i.e., if the 
difference between the atomic level spacing and the fre- 
quency of the light field is much larger than the coupling 



strength, the singularity at the energy shell is avoided, 
and the Frohlich transformation is expected to work well. 
It will result in an effective Hamiltonian which can also 
be obtained from the adiabatic elimination method Q 
and is equivalent to generic second-order perturbation 
theory. 

However, in quantum optics and atomic physics, most 
realistic systems involve time-dependent classical fields. 
Obviously, a time-independent Frohlich transformation 
cannot work for these cases. In this paper we use an 
effective Hamiltonian approach to eliminate certain in- 
termediate degrees of freedom (e.g., the photon degree 
of freedom, or the atomic operators concerning a certain 
atomic level) corresponding to the time-dependent terms 
for the general case. A time-dependent transformation 
is used to make the first-order terms zero and to keep 
the second-order terms. We will call this method time- 
dependent Frohlich transformation (TDFT). 

The paper is organized as follows: in Section [H] we will 
briefly introduce the TDFT to fix the notation. In Sec- 
tioii lllll and the following sections, we will use the TDFT 
approach to study the creation of two-atom entanglement 
after two atoms successively pass a cavity with a single 
mode field. As a result, we find that the entanglement 
depends on the atomic velocity (i.e., the transit time) 
and the initial distance between the two atoms. We de- 
termine the parameter regions for which the two atoms 
are maximally entangled, and discuss the limits of appli- 
cability of our method. 



II. TIME-DEPENDENT CANONICAL 
TRANSFORMATION 



We consider a general time-dependent quantum system 
with a Hamiltonian 



H = H U 



H\t), 



(3) 



where the unperturbed part H° is time- independent, 
and H x (t) is a time-dependent perturbation (|i? 1 (f)| <C 
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We perform the following time-dependent trans- 
formation: 

|^(t)) = e- s «Mf)>, (4) 

where is a quantum state whose time dependence is 

governed by the Hamiltonian H . The transformed state 
\ip'(t)) evolves according to the unitarily transformed 
Hamiltonian 

H' = er s He s + i (d t e~ s ) ■ e s 

oo 
n—l 

+ 4£^)-(££). (=) 

\n=0 / \n=0 / 

Keeping the terms up to second order, we obtain the 
following effective Hamiltonian: 

H' eS = H Q + (H 1 + [H°, S] - id t S) 

+ ^[{H l + [H°,S]-id t S),S] + \[H\S\. (6) 

If the operator S(t) in the time-dependent transforma- 
tion is chosen such as to make the first-order term of the 
effective Hamiltonian zero, that is 

H 1 + [H°,S] -id t S = 0, (7) 

the effective Hamiltonian takes the simple form 

H' eS = H°+~[H\S}. (8) 

The canonical transformation described by Eqs. ([7]) and 
([8]) is the so-called time- dependent Frohlich transforma- 
tion (TDFT). 

The form of the effective Hamiltonian after TDFT is 
similar to that in the time-independent case. Moreover, 
one also needs to have the formal solution for Eq. {7J to 
give an explicit expression for the effective Hamiltonian 
(0). We assume {\m) \m — 0, 1, 2, ...} to be a set of eigen- 
states for the time-independent zeroth-order Hamiltonian 
H°, and E m to be the eigenvalue for \m). The matrix 
elements of Eq. (J7J in this basis lead to 

H-mn + C^m — E n )S mn — id t S mn = 0, (9) 

where we have used the notation (dtS) mn = dtS mn since 
m) is time- independent. Thus, the solution for the trans- 
formation matrix elements S mn is given by @ 

= fe-^t-^H^dt'Win] (10) 

m,n ^° 

with E rnn := E m — E n for any m and n. 

We would like to remark that TDFT is equivalent 
to second-order time-dependent perturbation theory (as 



proved in Appendix YK\ . However, the effective Hamil- 
tonian obtained after applying TDFT, Eq. (jHJ), which 
contains only second-order interaction terms and no first- 
order interaction terms, is usually more convenient to be 
evaluated than the original Hamiltonian. For the time- 
independent weak-coupling atom-photon system (e.g., 
the systems given in Refs. [1, H3: or the electron-phonon 
systems in Refs. 0, Q), time- independent perturbation 
theory will give a second-order perturbation solution in- 
volving an infinite-dimensional Hilbert space. The equiv- 
alent conventional Frohlich approach will give a decou- 
pled effective Hamiltonian which can involve only the 
atomic part and can reduce to an analytically-solvable 
Hamiltonian in a finite-dimensional atomic Hilbert space. 
Similarly, in the time-dependent case, the interaction 
terms are time-dependent in the original Hamiltonian 
(]3]), and TDFT will give an effective (time-dependent) 
Hamiltonian, which may be decoupled from the photonic 
part and can be evaluated in the atomic Hilbert space. 
Therefore, applying TDFT can be much more efficient 
than using time-dependent perturbation theory. 

In the following sections, we will give an example for 
the power of TDFT, namely, a system of two atoms that 
cross a single-mode optical cavity. This is by no means 
the only example. We expect that many atom-light or 
electron-phonon models that have been studied by the 
conventional Frohlich approach (or other equivalent ap- 
proaches) to eliminate certain degrees of freedom, are 
suitable to be evaluated by TDFT if the interaction de- 
pends on time. To conclude: TDFT is powerful for atom- 
light or electron-phonon models whose time-dependent 
interaction is perturbative. 



III. TIME-DEPENDENT MODEL FOR TWO 
ATOMS PASSING A CAVITY SUCCESSIVELY 

We now use the TDFT method developed above to 
study a realistic physics problem in quantum optics. En- 
tanglement is a defining feature of quantum mechanics 
that has no classical counterpart. It is an interesting 
issue to entangle two atoms separated by a large dis- 
tance that have no direct interaction. Numerous propos- 
als have been made for entangling atoms trapped in a 
cavity/cavities 0, M El El M M HE E El El, 01 • 
Here, we will propose a scheme to entangle two identi- 
cal two-level atoms, which are not trapped in a cavity 
but pass the single-mode optical cavity sequentially in 
transverse direction, see Fig. [T] Our goal is to calculate 
the degree of entanglement between these two atoms after 
this process. The coupling of the atoms to the cavity field 
is position-dependent, and their motion therefore causes 
a time-dependent coupling. If the coupling energy is as- 
sumed to be much less than the detuning between the 
atomic transition frequency and the optical frequency, 
that is, the large-detuning condition is satisfied, we can 
use the TDFT to get an effective time-dependent Hamil- 
tonian, which only involves atom-atom interaction terms 
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and eliminates the optical field. 

A similar idea to entangle two atoms crossing a far-off- 
resonant single- mode cavity has been studied [12]. How- 
ever, there are significant differences between Ref. [l2j 
and our model: (i) In the proposal both two-level 
atoms enter (or leave) the cavity with the same velocity 
at the same time. In our model, the two atoms have dif- 
ferent initial positions, i.e., they enter the cavity with the 
same velocity at different times. In fact, we will study the 
degree of entanglement between the atoms as a function 
of the difference in initial position, (ii) In Ref. [lH, the 
coupling between the atoms and the cavity is assumed 
to be constant. Therefore, it is possible to obtain an 
effective Hamiltonian with a reduced atom-atom interac- 
tion (assuming large detuning) by eliminating the pho- 
tons (i.e., by means of the conventional Frohlich transfor- 
mation or other equivalent approaches). In reality, the 
coupling depends on the position the atom and has a 
Gaussian shape. Thus, the coupling is time-dependent 
when the atoms cross the cavity, and this is our motiva- 
tion to use TDFT to study the present model. There 
are many other works on the generation of entangle- 
ment between two atoms crossing a cavity. Reference 
[lij describes experiments that study entanglement be- 
tween two atoms after crossing a resonant-coupling cavity 
one by one. A two-qubit Grover quantum search algo- 
rithm is studied in Ref. [l5[ by looking at two atoms 
crossing a lar ge-d etuning cavity (this model is similar to 
that in Ref. [l2j ). The generation and purification of 
maximally entangled states of two A-type-atoms inside 
a large-detuning cavity have also been investigated [lg| . 
The important difference between all of the above models 
and our model is that the atom-photon coupling in the 
cavity is assumed to be constant in these models but is 
time-dependent in our model. 

As can be seen from Fig. [TJ the atoms are assumed to 
have the same constant velocity v along the z-direction 
with both the initial positions z® and z\ far away from 
the cavity. The Hamiltonian reads (H = 1) 

H or i = u)o(&i + ff) + wa'a 

+ gi{t)cr^a + g 2 (t)a 2 a + h.c, (11) 

where — \e)jj(g\ and crj = |e)jj( e l with |e) and \g) 
the atomic excited and ground states; gj(t) = g(z® + 
vt) sin(fc:z;) are the coupling constants which are assumed 
to be of Gaussian form, g(z) = g exp(— z 2 /d 2 ), where 
2d is the width of the cavity and go the vacuum Rabi 
frequency. Here we set sm(kx) = 1 since the atoms are 
assumed to pass the cavity in transverse direction at a 
maximum of the standing light wave. It is also assumed 
that the atomic velocity along the x- and y-direction is 
zero. We ignore the back-action of the cavity to the mo- 
mentum of the atoms since we assume the velocity to be 
large, (e.g., v ~ 10 m/s), such that the atomic decay time 
is long compared with the atomic cavity transit time. We 
assume the condition of large detuning is also fulfilled, 
|A| >• I <7o I ?i |<7i, 2 (*) I with the detuning A = ujq — oj, 




FIG. 1: (Color online) Two atoms crossing the optical cavity. 
The width of the cavity is 2d and determines the profile of 
the coupling constant g(z) = go exp(— z 2 /d 2 ). 

where ojo (uj) is the atomic transition frequency (optical 
frequency) . 

In the interaction picture with respect to 

H° ri =w(al + <4)+ua!a, (12) 

the Hamiltonian (fTTj) reads H = Ho + Hi (t) where 

Ho = A(of + erf), 

Hi = g 1 (t)cr+a + g 2 (t)a+a + h.c. (13) 

For large detuning, we can use a time-dependent Frohlich 
transformation defined by 

S(t) = xi(t)<rfa + x 2 (t)a 2 h a - h.c. (14) 

to eliminate the photon operators. The coefficients 
Xi, 2 (t) satisfy the following equations: 

g 3 {t) + Axj{t) - ixj(t) = 0, (j = 1, 2) (15) 

according to Eq. ([7]) . The explicit solution for the coeffi- 
cients xi^(t) is given by [20j 

xj(t) = -i fgtfyW-Vdt', (j = 1,2). (16) 
Jo 

After the Frohlich transformation, the effective Hamil- 
tonian reads (up to a constant term) 

H' eS {t)=H + ±[H 1 (t),S(t)] 

= A lC r? + A 2 <7| + (7+ + / V^ffJ, (17) 

where 

A, = A - (g* Xj + 9j x*)(l + 2n p ) (18) 

for j = 1,2, rip = (a^a) is the mean photon number in 
the cavity, and 

f(t) = -\{glx 2 +g 2 x\) (19) 
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is the effective coupling between the two atoms. The 
(first-order) weak interaction terms between the photons 
and atoms have been eliminated, and only the induced 
(second-order) atom-atom interaction terms remain. The 
Hamiltonian (|17p . which contains no photon operators 
(except the time-independent expectation value of the 
photon number operator), is much simpler than the orig- 
inal one (fTTj) . It can be diagonalized in the atomic basis 

{\9192), \g1e2), \eig 2 ), |eie 2 » = {\gg), \ge), \eg), \ee)}. 

(20) 

So far, we do not need any condition except the large 
detuning to get the simple effective Hamiltonian (|17|l . In 
the following section, we will further simplify the effec- 
tive Hamiltonian (|17[) using typical parameters for the 
cavity-QED system. As we will see below, the adiabatic 
condition is satisfied since the time-dependent couplings 
9i,2(t) are slowly varying (that is, \gj(t)/gj(t)\ <C |A| for 
j = 1,2). Of course, if the adiabatic condition is satisfied, 
the result obtained by TDFT in the next section can also 
be obtained by using adiabatic elimination. However, the 
elimination approach cannot be expected to describe our 
model if the adiabatic condition is not fulfilled. 



IV. DYNAMICAL GENERATION OF 
TWO-ATOM ENTANGLEMENT 

In the following, we will simplify the Hamiltonian (JTTJ) 
further by taking into account typical parameters for the 
cavity-QED system. 

First, we will rewrite the coefficients ([Till) as follows: 



Xj(t) 



-le 



-iAt 



9 j(t')e iAt 'dt' 



g 3 (t')e^ dt 



(21) 



Here, we have used the explicit Gaussian form of the 
coupling constants. We have also used 



(z° + u%(tV dt' -> d / 5j(^)e dt' (22) 



since z° + vt ~ d in the effective integration range and 
the change of gj (t) (also change of z° + vt) is much slower 
than that of e iAt : 

F{g 3 {t)) ~ F(*° + vt) ~ I F(e-^) ~ A, (23) 

where F(x) := dtx/x denotes the change of x. For typical 
system parameters [2l|, [22|, [23[ (wavelength of the atomic 
transition A = 1000 nm, detuning A = 10 4 MHz, go — 
100 MHz, cavity width (in the z-direction) d = 30A, and 
atomic velocity of the order of v = 10 m/s), (v/d) ~ 
10 _1 MHz < |A| = 10 4 MHz, that is, the atomic cavity 



transit time t trans i t = 2d/v is much larger than the time 
1/A according to the detuning. Actually, the fact that 
v/d < |A| (that is, \gj{t)/gj{t)\ < |A| for j = 1, 2) means 
the adiabatic condition is satisfied. 

Using these parameter values, Xj(t) in Eq. (f2~Tj) can be 
further simplified as 



/(*) 



<?;(*) 



lQ- 5 ixj(t) 



The coefficients in Eq. (JTTJ) thus take the following 
form: 



A, 



2{l + 2n p )g](t) 



-J - " -r A 

and the effective coupling is 



/(*) 



A 



(24) 



(25) 



We will now further assume that the shift terms in 
Aj can be ignored since 2g|(i)/|A| <C |A|, i.e., we will 
replace the coefficients Aj by A. Hence, 



H' eS (t) = Aa{ + Aa z 2 + fa^a+ + fa+a^ 



(26) 



We now apply the transformation U(t) = exp{(Acrf 4 
Aa%)t} to the Hamiltonian H' eS (t) defined in Eq. JH 
and obtain 



H' c ' s (t) = fa^a+ + f f j+a-. 



(27) 



The time evolution of a general state governed by H' c ' s (t) 
is given as 

m)) = C gg {t) \gg)+C ge (t) \ge)+C eg (t) \eg)+C ee (t) \ee) , 

(28) 

where C gg (t) = C gg (0), C ee (t) = C ee (0), and 

C ge (t) = C ge (0)cos9(t) -iC eg (0)sm0(t), (29) 
C eg {t) = -iC ge {0)sin9{t) +C eg {0)cos9{t), (30) 

with 6{t) = jlf{t')dt' ~ J^fWdtf. Hence, the sys- 
tem described by Eq. (|2T|) is characterized by a closed 
subspace {\ge),\eg)}. 

In what follows, we will investigate what degree of en- 
tanglement can be obtained after both atoms pass the 
cavity in transverse direction. We assume that the two 
atoms are prepared in the initial state |"0(O)} = \ge). The 
time evolution generated by Eq. (|27p leads to 



|VC*)> = cos 9(t)\ge) - i sin 6 (t)\eg). 



(31) 



After both atoms have passed through the cavity and are 
far outside, that is t — > 4- 00, one has 



(z°) 2 



2d 2 



(32) 



where z° = z\ — z° denotes the difference between the 
initial atomic positions. 
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FIG. 2: (Color online) Entanglement entropy E as a function 
of atomic velocity v (in units of gfid/A) and initial atomic 
position difference z° (in units of d). 

In general, the state shown in Eq. (j3"Tj) corresponds to 
an entangled state of the two atoms. In the following, 
we will study the degree of entanglement as a function 
of z° and the atomic velocity v using the entanglement 
entropy that is defined as 

£(|^))) = -THpilog 2 pi(t)). (33) 

Here, \ip(t)) is a pure state, and pi{t) = Tr 2 (|^(i)) (V'(OI) 
is the reduced density matrix of the first atom. Evalu- 
ating this expression for the state shown in Eq. (|31[) . we 
obtain 

£(|V>(+oo))) = - cos 2 [0(+oo)] log 2 cos 2 [0(+oo)] 

- sin 2 [9(+oo)\ log 2 sin 2 [0(+oo)] . (34) 

A maximally entangled state for the atoms occurs for 
8(+oo) = (2n + 1)tt/4 for any integer n. Figure [2] shows 
the entanglement entropy E(\ip(+oo))) as a function of 
initial atomic position difference z° (in units of d) and 
atomic velocity v (in units of g^d/A, which is 30 m/s 
for the system parameters discussed after Eq. (p23"|) ). The 
possible values for z° and v which make the entanglement 
maximal (that is, E(\ip(+oo))) = 1) are shown in Fig. [3] 

It is interesting that there may be maximally entan- 
gled states for the two atoms even if the first atom has 
left the cavity before the second atom begins to enter it. 
For example, when the velocity v is smaller than about 
0.23g 2 <i/A (corresponding to about 7 m/s for the sys- 
tem parameters discussed after Eq. (|23|) ), it is possible to 
obtain a maximally entangled state although the initial 
atomic position difference z° can be larger than 2d (i.e., 
the approximate transverse width of cavity), see Figj3] 

In the above calculation, we have ignored the effect of 
(photonic and atomic) decay and photonic back-action to 
the atoms for the following reasons: (i) we assume that 
the velocity is of the order of 10 m/s. Higher velocities 
will lead to a reduced degree of entanglement; lower ve- 
locities will lead to optical back-action to the atoms and 
photonic/ atomic decay, (ii) we have assumed large de- 
tuning to reduce the atomic decay, (iii) the mean number 
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FIG. 3: (Color online) Relation between atomic velocity v (in 
units of Qgd/A) and initial atomic position difference z° (in 
units of d) along the lines of maximal entanglement E(+oo) = 
1, or 0(+oo) = (2n + l)vr/4, for n = 0, 1, 6. 

of photons n p is assumed to be small (e.g. 10 -4 ) in order 
to reduce the effect due to photonic decay. For exam- 
ple, the atomic cavity transit time (which is close to the 
interaction time) is less than the effective photonic deco- 
herence time: 2d/v ~ 10~ 5 s < \/(n p n) ~ 10~ 4 s (where 
k ~ 10 8 Hz is the typical cavity one-photon damping rate 
|22|). Also, we have chosen Ai ;2 (i) — > A which amounts 
to ignoring the effect of the shift terms in Eq. (|24jl . The 
explicit form of Ai i2 (t) shows that if the initial atomic 
position difference z° is close to 0, then Ai(t) — * A 2 (t), 
and the effect of the shift terms can be neglected. 

Finally, we would like to remark that the time- 
dependent Frohlich transformation is valid under the con- 
dition of weak coupling (which is equivalent to the con- 
dition of large detuning here): H\{t) H n (g A). 
The adiabatic condition that the atomic cavity transit 
time ttransit = 2d/v is much larger than 1/A, is only 
used to simplify the effective coupling in Eq. (fT7|) for 
the present realistic atoms-cavity system. Actually, this 
adiabatic condition is independent of the above large de- 
tuning condition for the time-dependent Frohlich trans- 
formation. The above argument proves that, the TDFT 
approach is valid independent of whether the adiabatic 
condition is fulfilled or not. 



V. CONCLUSION 

Using the time-dependent Frohlich transformation, we 
have calculated the degree of atomic entanglement be- 
tween two atoms passing an optical cavity sequentially. 
The Frohlich transformation eliminates the photonic op- 
erators and induces an effective atom-atom interaction. 
We have determined the velocities and initial atomic po- 
sition differences for which the entanglement is maximal, 
and we have shown that there may be maximally entan- 
gled states for the two atoms even if the first atom has 
left the cavity before the second atom begins to enter it. 

A number of time-dependent canonical transformation 
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methods have been proposed 0, [24| to discuss time- 
dependent problems. The difference of our approach to a 
general time-dependent canonical transformation is that 
our TDFT requires the generator S(t) in the canoni- 
cal transformation to satisfy Eq. ||7J) in order to elimi- 
nate some intermediate degree of freedom for a general 
system. As an example, we have studied a system of 
two atoms that pass a cavity, which can be described by 
an effective atom-atom interaction Hamiltonian by using 
the TDFT to eliminate the photon degrees of freedom. 
In fact, the TDFT can be used to treat certain time- 
dependent atom- light systems by eliminating some of the 
atomic degree of freedom (e.g., one atomic level) like the 
Frohlich transformation for the time-independent cases. 
The TDFT presented here works well not only for opti- 
cal and atomic systems, but also in the field of condensed 
matter physics (e.g., time-dependent electron-phonon in- 
teractions) or other systems with weak (first-order) time- 
dependent interactions. 

This work was supported by the European Union under 
contract IST-3-015708-IP EuroSQIP, by the Swiss NSF, 
and the NCCR Nanoscience. It was also supported by 
the NSFC with grant Nos. 90203018, 10574133, 10474104 
and 60433050, and NFRPC with No. 2005CB724508. 



APPENDIX A: EQUIVALENCE TO 
TIME-DEPENDENT SECOND-ORDER 
PERTURBATION THEORY 

To check the range of validity of the above effective 
Hamiltonian obtained by TDFT, we compare with the 
results of standard time-dependent perturbation theory 
up to of second order. 

To apply time-dependent second-order perturbation 
theory for the Hamiltonian H = H° + H 1 ^) given in 
Eq. ([3]), we assume { | m) \m = 0, 1, 2...} to be a complete 
set of eigenstates for the time-independent zeroth-order 
Hamiltonian H°, with eigenvalues E m . A state \ip(t)) 
satisfies the Schrodinger equation 

id t m)) = [H°+H i (t)}\m)- (ai) 

We now do perturbation theory in H 1 (t) by assuming 
\m) = J2\^ l) ®)= E CW(t)e-^-*|n). (A2) 



1=0 



=0,1=0 



Here, \if)( l \t)) is the Z-order contribution to the pertur- 
bation expansion. Replacing \ip(t)) in Eq. (|Aip by the 
perturbation expansion and comparing the coefficients 
order by order, we obtain the following equation: 

*Ci +1 \t) =J2c^(ty E - t HUt)- (A3) 



If the initial state is \k), the zeroth-order solution for 
the coefficients C n (t) is 



cl 0) (*) 



(A4) 



According to Eq. (IA3|) the first-order solution is 

ft 



-lie 
'o 



iE m kt' ttI UUJ 



(A5) 



The second-order solution is 
rt 



(A6) 

We will now consider the time evolution of the state 
predicted by the TDFT method. It follows from the 
Schrodinger equation 

idtW{t)) = {H°+ l -[H\S])W(t)) (A7) 

governed by the effective Hamiltonian H' eS in ([8]) after 
applying the TDFT method to the original Hamiltonian 
H = H° + H x {t). 

The connection to the original picture before the 
TDFT, i.e., the state |^>(i)), whose time dependence is 
governed by H, reads 



|0(t))^(l + 5 + -5 2 + ...)|^(i)). 



(A8) 



Like in Eq. (|A2[) . we can write down the perturbation 
expansion for if)'(t): 

|V/(i))=El^ (0 (i)>= E C'^(t)e- iE ^\n). (A9) 



1=0 



=0,1=0 



Starting with the initial state \k), and keeping terms 
up to the zeroth-order in Eq. (|A8j) . it is obvious that 



'(0), 



and 



C<®(t) = C'W(t) = S mh . 



(A10) 



(All) 



Since the perturbation term in H' eS is of second order, the 
first-order correction for \ip'(t)) vanishes: \if)'^(t)) = 0. 
Then from Eq. (|A8|) . the first-order correction for the 
state \ip(t)} is 



(A12) 



Correspondingly, the first-order correction coefficients 
are 



<4 x >(t) 



-lie 
lo 



iE mk t' rjl 



H x mk (t')dt'. 



(A13) 



The results of Eq. (|A13[) are equivalent to the stan- 
dard time-dependent perturbation theory as shown in Eq. 

(EH). 

The second-order correction for \ip(t)) is 



(A14) 
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where \ip ,l(2 \t)) is the correction according to [H 1 ,S]/2 
in H' cff © with the exact form being given as (according 
to Eq. JA5|) 

^\t))=Y,C'i 2) {t)e- iB ^\m) (A15) 

m 

with 

C'i 2) (t) = -\ j\ iE ^'[H\t%S{t')] mk &t>. 

Then the coefficients for the second-order corrections of 
the state |^ (2) (*)) are given as 

C${t)=C'W{t) + \Y. S ™nSnk. 



By replacing the operator S(t) by that in Eq. (|T0|) in the 
above equation, we obtain 

n J ° 

x /' At"e iE ^"H l nk {t'% (A16) 
Jo 

which agrees with the expression given in Eq. (|A6[) by 
time-dependent second-order perturbation theory. 

In this appendix, by reformulating the well-known per- 
turbation solution to the time-dependent Schrodinger 
equation in our notations, we have shown explicitly 
that the second-order time-dependent perturbation the- 
ory agrees with TDFT. 
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